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Abstract 


A mathematical model is presented for a lithium ion cell under galvanostatic discharge. The non-homogeneous material balance equation for 
the solution phase of a lithium cell was solved analytically using Green’s function. We determined the elementary solutions, the norms and the 
eigen values of the problems for galvanostatic boundary conditions and systematically tabulate the resulting expressions. The solution has been 


obtained in a closed form expression. 
© 2007 Elsevier B.V. All rights reserved. 


Keywords: Green’s function; Non-homogeneous equation; Concentration; Eigen values; Eigen functions; Utilization; Boundary value problem 


1. Introduction 


The use of mathematical modeling in the design of battery 
systems has a long history. Macroscopic models of the current 
and potential distribution in porous battery electrodes were first 
developed in the late 1950s [1,2]. Major strides towards under- 
standing the behaviour of porous electrodes were made in the 
early sixties with the development of porous electrode theory 
[3] which generalized earlier modeling efforts into a macro- 
homogeneous modeling framework that is still used in most 
present day models. In the early seventies, the scope of battery 
modeling increased from treating a single electrode to modeling 
the full-cell sandwich including two electrodes and the separa- 
tor by numerical computation [4,5]. This modeling approach 
allowed the treatment of complicated interactions between a 
wide variety of phenomena that had previously been studied 
separately under limiting conditions. 

The earliest model for lithium-ion processes in the compos- 
ite insertion electrode was developed in the 1980s [6-8]. The 
model covers only a single porous electrode. It does not have 
the advantages of a full-cell sandwich model for the treatment 
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of complex interacting phenomena between the cell layers and 
utility for design purposes. Doyle et al. [9,10] have developed a 
very general one-dimensional model applicable to almost any of 
the existing Li/Li* systems. Their model uses the porous electro 
theory developed by Newman [5] to describe the potential varia- 
tions in the solid and solution phases. The material balance in the 
solution phase is described using the concentrated solution the- 
ory and in the solid phase using a Fickian diffusion equation in 
spherical coordinates. Confining only to thermal conditions, the 
authors validated their model by demonstrating good agreement 
with experimental data. Later, this model has been extended to 
include an energy balance in order to predict cell temperature 
[11-13]. The problems encountered with Li-ion systems like 
poor rate capability, thermal runaway, occurrence of undesired 
side reactions and capacity fade have been addressed in 2002 
using the new model [14]. 

Numerical techniques such as the finite element method and 
finite difference method allow results of a high accuracy but 
they demand extensive computing time. Analytical solutions 
are the best if available due to the solutions being continu- 
ous in the independent variables and show explicitly how the 
parameters of the system are involved. The solutions give much 
insight into a system, which is one of the primary objectives 
of modeling. Unfortunately, analytical solutions are not avail- 
able for most of cases of interest. It is specific to the system, 
geometry and boundary conditions. Usually, analytical solu- 
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Nomenclature 

a specific interfacial area (m? m~?) 

c dimensionless concentration 

Ci concentration of species i (mol m7?) 

Diett effective coefficient diffusion of species i 
(m? s7!) 

F Faraday’s constant (996,487 C cq7 !) 

I superficial current density (A m~7) 

jn pore-wall flux of lithium ion (mol m~? s7!) 

J dimensionless flux 

L thickness of the cell (m) 

Le thickness of the cathode (m) 

Ls thickness of the separator (m) 

Ni molar flux of species i (mol m~? s7!) 

r dimensionless thickness of the cell/separator 

t time (s) 

ti transference number of species i 

x distance from the separator/positive electrode 
boundary (m) 

y dimensionless distance along x-axis 

Zi charge per mole 


Greek letters 


E porosity of electrode 

Àn eigen value 

T dimensionless time 

W dimensionless concentration (for homogeneous 
equation) 

Subscripts 

I We Die sca 

n 0,1,2,... 

eff effective 

+ positive electrode 

Li lithium 

c cathode 

s separator or solid phase 


tions are restricted to linear equations with linear boundary 
conditions. In 1979, Atlung et al. [6] solved the diffusion equa- 
tion for solid solution cathodes with zero initial concentration 
using Laplace transformation technique. In 1997, Doyle and 
Newman [10] solved a similar model for the discharge lithium 
batteries under steady-state condition using the method of sep- 
aration of variables. In 2001, Subramaniam and White [15] 
used an extended method of separation of variables for different 
electrode geometries with galvanostatic boundary conditions at 
zero initial concentration. This was followed by the work of 
Hashim et al. [16,17] who solved the material balance equa- 
tion in lithium-ion batteries with non-zero initial concentration 
by using Laplace transform. In 2004, Johan and Arof [18,19] 
presented analytical solutions for lithium intercalation in dif- 
ferent cathode geometries using integral transform technique. 
Recently, Liu [20] presented an analytical result of the intercala- 


tion process of Li/Li* into a spherical particle of graphite/cobalt 
oxide immersed in a conductive electrolyte using finite integral 
transform method. In his work, Liu applied pseudo-steady-state 
condition in order to guarantee a uniform convergence of the 
solution. 

The objective of this paper is to demonstrate a method 
of Green’s function in solving material balance equation in 
solution phase of electrochemical system. The Green’s func- 
tion approach is a powerful technique for finding solutions 
to linear, homogeneous and non-homogeneous partial differ- 
ential equations [21-27]. To our best knowledge, Green’s 
function approach has not been used for the lithium battery 
systems. 


2. Specific capacity of lithium ion cell 


In this work, a model has been developed for a galvanos- 
tatic discharge behaviour of a Li/LiMn20; as shown in Fig. 1. 
During the discharge process, lithium ions de-intercalate from 
the anode diffuse through the separator and intercalate into the 
cathode. For this work, the cathode is closely packed electrode 
employing a LiMn20;4 active material. For one-dimensional 
transport of lithium ions in a porous composite cathode with 
active materials, the material balance equation for species i an be 
written as 


d(Ec;) 
ot 

where € is the porosity and is considered to be constant, c; the 
concentration of the species i, N; the species flux, a the species 
interfacial area, jn the pore wall flux across the electrolyte active 
material interface and t; is the transference number of species i. 
From the concentrated solution theory [5], the series flux is 

tT 

Zi 


==VNi täh — ti) (1) 


Ni = —£ Di effV ci + 


(2) 


where J is the electrical current density and z;F is the charge 
per mole. The term Dj efr is the effective diffusion coefficient of 


anode 
cathode ——} 


separator 


Fig. 1. Schematic representation of lithium-ion battery. 
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species i defined by Dj ett = Dye?!” [17], where D; is the diffusion 
coefficient of species i in the electrolyte. Substituting Eq. (2) into 
Eq. (1) gives 

Ochi It; 


E 
ot 


L 


= —V(—e Di eff V ci) — V ( ) + ajn( — ti) (3) 
Since i is a lithium ion species, then t; = ti = t+ =0.2 [10]. We 
can further simplify Eq. (3) as 


Tiia ajn (0.8) 
3x? E 

The distribution of reaction rate in the porous electrode is 
a complex process. It reflects the solution-phase concentration 
variations as well as the trade-off between ohmic and kinetic 
resistances. When using an insertion material having an open 
circuit potential that depends strongly on the state of charge of 
the system, a uniform current distribution results [28]. Other than 
that, a uniform current distribution can be expected when kinetic 
resistances dominate ohmic resistances [5]. According to [10], 
for a uniform current distribution, it is possible to assume that 
the pore wall flux across the electrolyte—intercalation material 
interface is given by its average value at any point in the porous 
electrode. Hence 


I 
aFLe 


(4) 


Jn = (5) 


where Le is the thickness of the cathode and 7 is the discharge 
current. Introducing the dimensionless parameters 


CLi x L Dy jt 


C=; y= —, f= T= (6) 
CLi,O Ls Ls L2 

Substituting Eq. (6) into Eq. (4) gives 
ð 3? 
Zap" 45 p0 1<y<ltr (7a) 
at dy 

Subjected to the boundary and initial conditions 
ð 
£0 230, 12921 42 (7b) 
dy 
dc Jr 
dy el?’ t>0, y=1 (7c) 
c(y, 0) = 1 (7d) 
where dimensionless flux is 

0.8)L2j 
J= a(0.8)LS J (8) 
EDCLi,O 


We solved Eq. (7) by using Green function. The homogeneous 
version of Eqs. (7a)—(7d) are: 
ey 1 dw 


— = t> 0, 


ye ele ae I<y<l+r (9a) 


subject to the boundary and initial conditions 


ve 


=0, t>, y=l+r (9b) 
dy 


Map. t>0, y=1 (9c) 
dy 
y, 0) = FY) =1 (9d) 


The solution for the homogeneous Eqs. (9a)—(9d) is 


Wy, 1) = YAn X On, ype" (10) 


n=0 
where the unknown constant given as 


1 
N(An) 


n= 


l+r 
f XOn VFO dy’ a1) 


and the norm is 


l+r 2 
Nn) = f [X@n, yT dy (12) 


Hence, Eq. (10) becomes 
co Aer 


eo An 1+r 
HW. = Xn”) | Xn, Y)FOAY (13) 
n=0 


Rearrange Eq. (13) 


[tr | œ eit 
Wy, T) = f È Xon 0m | Fo’) dy’ 


NOn) 

(14) 
where eigen functions given as 
X(An, y) = cos AS) (15a) 
Xo, y) = 1 (15b) 

and the norms are 

NOn) = 5 (16a) 
No) =r (16b) 


From Eq. (14), we an determine for the homogeneous Green’s 
function as 


—_— = (ltr = 2 
4 ZS M ieee te Ant 

Oto. ry¥.0) = 7 + 7D eos ( AT e 
n= 


— (£ N ) a7) 


el/4 


Hence, non-homogeneous Green’s function given as: 


el/4 


a (=£ cha z) (18) 


el /4 


i 2 An +r—y a 
GO, tly, “)=-+-— | ea 
(y, tly’, t’) 7+ Fem ( e 
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The general solution for Eq. (7) is 


l+r 
c(y, t) = f G(y, tly’, OFO" 


r l+r 
a f ae f G(y, tly’, dy’ 
0 il 


a 1+ 
+l yr | ar f GO, TY’, Dyan a9) 
0 1 


Substituting Green’s functions Eqs. (17) and (18) into Eq. 
(19) and perform the integration, finally we obtain 


Ji 1 
c(y, T) = 1 + Jrre? + Jre? + > (> r +r 3 ) 
r 


Àn(l — y 2 
x COS SS) evant (20) 
el/4 


where àn =nz/r. The term in the summation is only significant 
forn=1. 


3. Results and discussion 


The profile of the concentration of lithium ion during a gal- 
vanostatic discharge in the solution phase of the cathode is 
plotted in Fig. 2 using Eq. (20). Fig. 2 shows that the lithium 
ion concentration is minimum at the back of the cathode, that is 
at y=1+r. This is reasonable because during full charge, it is 
expected that the lithium ion concentration in the cathode will 
be almost zero. During discharge, the lithium ion will fill up the 
back portion of the cathode last. The result displayed in Fig. 1 
shows a similar trend to that obtained by Doyle and Newman 
[10], Hashim et al. [16,17] and Subramaniam and White [15] 
as shown in Fig. 3. The intersection point of Fig. 2 occurs near 
y=3 while in the graph displayed in reference [10] show inter- 
section at y=2 and in Ref. [15] show intersection at y= 1.5. 
This is due to the effect of taking the boundary condition at the 
separator/cathode interface as Jr/e!/? while Doyle and Newman 
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Fig. 2. Concentration profile in solution phase across cathode during galvanos- 
tatic discharge as a function of applied current at a particular time (T = 1). 
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Fig. 3. Concentration profiles in solution phase across cathode compared with 
published results. 


have assumed a steady-state situation when plotting this graph. 
The result displayed in Refs. [16,17] show a similar intersection 
point (y= 3) with ours. 


4. Conclusion 


This paper presents the analytical solution for the concentra- 
tion profile of lithium ion in the solution phase of the positive 
electrode of a lithium ion cell. The model is justified for solu- 
tion phase diffusion limitations in the discharge systems. The 
differential and algebraic equations that describe the concen- 
tration, potential and current density in the solution phase of 
the porous electrode are in general nonlinear and coupled [9]. 
To find an analytical solution we must first decouple these 
equations by assuming a form for the reaction rate distribu- 
tion through the porous electrode. We made an assumption that 
jn is a constant at any point in the porous electrode. The con- 
centration profiles will be determined using the boundary and 
initial conditions. The method Green’s function presented here 
is useful for solving boundary value problems that include the 
flux boundary conditions. This method yields an unambiguous, 
straightforward way to obtain analytical solutions for lithium 
ion diffusion. In solving for the concentration profile, the initial 
concentration is taken to be cg instead of a zero initial condi- 
tion as used by both Carslaw and Jaeger [27] and Atlung et 
al. [6]. The method applied can be easily extended to spheri- 
cal and cylindrical coordinates for concentration profile in solid 
phase. The method appears to be more general, systematic 
and straightforward and should be useful for solving the sim- 
ilar problems. However, a more detailed mathematical model 
is generally required for the final optimization of the battery 
design. 
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